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Abstract 

A new method based on modified optical Bloch equations is proposed for 
studying resonant tunneling in semiconductor heterostructures. The method 
manifestly takes account of electrons' statistics, which enables to investigate 
the influence of the Pauli principle on resonant tunneling in presence of inelas- 
tic scattering. Being applied to evaluation of the resonant current in semicon- 
ductor heterostructures, our approach predicts considerable deviations from 
the one-electron picture. 
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The impressive progress in microfabrication technology has considerably enhanced the 
interest to quantum transport through multilayered semiconductor heterostructures A 
good number of theoretical approaches have been developed to describe resonant tunneling 
in these devices. However, the problem has been mostly treated within the one-electron 
picture. The existing many-electron approaches (see, e. g., are very complicated, 

owing to which the influence of many-body effects in resonant tunneling problems remains 
insufficiently elucidated yet. One of them is the effect of the Pauli exclusion principle, which 
is expected to be especially important in presence of inelastic scattering 

We propose here a novel method of handling resonant-tunneling problems. This method 
is based on modified master equations suited for problems, where "classical" (incoherent) 
and "quantum" (coherent) phenomena are in interplay. Being technically simpler than the 
other approaches, it has nevertheless a many-body nature and takes account of the Pauli 
principle from the very beginning. 

Let us consider a mesoscopic "device" consisting of coupled quantum wells. The device 
is connected with two separate electron reservoirs. We assume that the density of states in 
the reservoirs is very high (continuum). However, each of the quantum wells in the device 
contains only a small number of discrete energy levels. For the sake of simplicity we consider 
the reservoirs at zero temperature, where the carriers (electrons) occupy all the states in the 
left reservoir up to the Fermi level e'P and in the right reservoir up to the level e'P < Ep\ 
As a result, a current is fiowing through the device from the left reservoir to the right one. 
The device is also coupled to an additional phonon reservoir (also at zero temperature), 
which enables inelastic transitions between the discrete levels in the quantum wells. We 
assume that the system is described by a tunneling Hamiltonian, i.e. only the couplings 
between the adjacent wells are taken into account. 

The general idea of our approach follows from the observation that the electron transport 
in the system can be fully determined in terms of the density submatrix of the device 
solely, without resorting to the complicated density matrix of the entire system including 
the reservoirs. This density submatrix contains full information about currents and particle 
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accumulation, and the only problem is to determine its time evolution. This problem can 
be solved by using, with some modifications, the technique of NMR theory ^ or the Bloch- 
Maxwell equations 0. 

Let us assume that the reservoirs' states remain unchanged during the process. We 
designate the density submatrix of the device as (Jab{t) in terms of all possible electrons' 
states (I a), |c), . . .) inside the device. For instance, \a) denotes the state where there are 
no carriers inside the device; \h) - one electron occupies the upper level in the first well, and 
so on. In the spirit of Bloch-Maxwell equations we write the equations for diagonal and 
off-diagonal matrix elements of the density submatrix as 

C^aa = -i\H,a]aa - aa ^ ^ a-*n + ^ (T c^a , (la) 

n(^a) c{j^a) 



O-ab = -i['H,(r]ab - ^CTab ( ^a^n + ^b^n | , (lb) 

\n(^a) n(^6) 



where Ti. is the matrix of the "internal" (without reservoirs) Hamiltonian, and Ta^n is 
the probability per unit time for the system to make a transition from the state \a) to 
the state \n) of the device due to the tunneling to (or from) the reservoirs, or due to 
interaction with the phonon bath. The commutator in these equations describes the coherent 
("quantum") motion of carriers inside the device, and the remaining terms describe the 
incoherent ("classical") motion of the carriers. By neglecting the commutator one obtains 
the "classical" kinetic master equations. 

Now we can use our basic assumption that the system is described by a tunneling Hamil- 
tonian, which takes into account only the couplings between the nearest neighbors. Then 
we can rewrite Eqs. (1) explicitly as 

O'aa = i X! ^abicTab " O'fea) " aa ^ ^a^n + ^ CTcc^c^a , (2a) 
b{j^a) n{y^a) c{^a) 

dab = i{,Eb — Ea)(Tab + i \ X! O^ab'^b'b— X! ^aa'^^a'i 

\b'{^b) a'i^a) 



Ta-.n + Y ^b-.n ] CTab , (2b) 
yn{j^a) n(^b) 



where Qij denote the couphngs between the levels in adjacent wells, and a^a = a*^- 

By solving Eqs. (0), one can find the diagonal elements of the density submatrix, which 
determine the charge density of the corresponding quantum states of the device. Then the 
current fiowing into the right reservoir is 

^W = E^^^r«, (3) 

i 

where the sum is taken over only those states \i) that contain electrons in the well adjacent 
to the right reservoir; T^^ is the partial width of the state \i) due to tunneling to the right 
reservoir |Q. The stationary current / is obtained from Eq. as / = J{t oo) and does 
not depend on the initial state of the device [P,p!0|] . 

Now we give some examples of application of our method to resonant tunneling problems. 
Although the method can be applied to many different problems, we concentrate in this 
letter on the Pauli exclusion effects in resonant tunneling. In all the examples we neglect 
transitions to motion parallel to the layers, by assuming that the scattering on impurities and 
layers irregularities is small. This makes the problem essentially one-dimensional. For the 
sake of simplicity we also neglect the electron spin effects, assuming only that two electrons 
cannot occupy the same quantum state. 

Let us start with the well-known case of penetration through a double-barrier heterostruc- 
ture with one discrete level Ei inside the well. Fig. 1(a). In this case the device has only 
two possible states, namely: \a) - the level Ei is empty, and \b) - the level Ei is occupied. 
There are two possible interstate transitions: \a) \b) via F^^ (an electron tunnels from 
the left reservoir), and \b) \a) via F/j (an electron tunnels to the right reservoir). Hence 
Eqs. (|2^) give 

0"aa = -FiCTaa + ^RC^bb , (4a) 
(^bb = ^L<^aa - FijO-ftfe , (4b) 

where Ti pi are the partial widths of the level Ei due to tunneling to the left and the 
right reservoirs. Solving Eqs. (§) with the initial condition either cTaa(O) = 1, Cbb = 0, or 
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o"aa(0) = 0, abb = I, and using Eq. we obtain for the resonant current / the standard 
result 



Notice that the coherent (quantum) effects in the resonant current are absent here. As it 
follows from Eqs. such effects can appear only due to coupling between discrete levels 
in adjacent quantum wells, whereas in the considered example there is only one quantum 
well. Hence, the transport through this system is described by classical master equations. It 
thus might be not surprising that "coherent" and "sequential" approaches to the resonant 
tunneling produce the same answer in this case, 0]. 

Next we consider the resonant tunneling through double-well potential structure shown 
in Fig. lb. As in the previous case, there are no inelastic transitions due to interaction with 
the phonon reservoir, since each of the wells contains only one level. However, the coherent 
tunneling between the wells can take place in this structure. Let us enumerate all possible 
electron states in the device. There are just four of them: \a) - the levels Ei^2 are empty; 
1 6) - the level Ei is occupied; |c) - the level E2 is occupied; \d) - both the levels Ei and 
E2 are occupied. Taking account of all possible transition between these states, one obtains 
from Eqs. 



<^aa = -'^L'^aa + ^RO-cc , (6a) 

abb = ^L(^aa + ^R(^dd + i^{(^bc " C^cfe) , (6b) 

O-cc = -^RCTcc - ^LCTcc + i^{cFcb " Cfbc) , (6c) 

^dd = -^RC^dd + TlcTcc , (6d) 

(Jbc = i{E2- Ei)abc + i^{(Jbb - c^cc) - ^ (Fl + ^r) (Jbc ■ (6e) 

The total current J, according to Eq. (|^), is / = {acc + o'dd)^R, where an = aii{t oo). 
Eqs. (^ can be easily solved, so that we finally get for the resonant current 



where e = E2 — E1. This result coincides with the one found in the framework of one-electron 
approach [ p!0| , p]2[ . At first sight, such a coincidence looks surprising. Indeed, our treatment 



explicitly excludes the states when two or more electrons occupy the same level. On the 
other hand, these states are implicitly present within the one-particle picture. Nevertheless, 
a more detailed analysis shows that this coincidence is not accidental. One can demonstrate 
that in the absence of inelastic scattering the Pauli forbidden configurations in the one- 



electron picture do not contribute to the resonant current [I^. However, in the presence of 
inelastic scattering, the infiuence of Pauli forbidden terms is very essential. 

Let us consider two other examples of the resonant tunneling, now in presence of inelastic 
scattering. The first case is shown in Fig. 2a. It corresponds to the double-barrier structure 
discussed above, but now the quantum well contains two levels, the upper one Eq and the 
lower one Ei. An electron which tunnels from the left reservoir into the upper level Eq can 
either relax inelastically into the lower state with the transition rate Fj^, and then tunnel 
out into the right reservoir, or tunnel out directly into the right reservoir. There are four 
possible electron states of the device: \a) - the levels Eq i are empty; \b) - the upper level, Eq, 
is occupied; |c) - the lower level, is occupied; \d) - both levels Eq and Ei are occupied. 
The possible interstate transitions are: \a) \h) via F/^, which is the tunneling rate from 
the left reservoir; \b) \a) via F^jo, which is the tunneling rate to the right reservoir from 
the upper level; \b) \c) via Fj„; |c) \a) via F/ji which is the tunneling rate to the 
right reservoir from the lower level; |c) — ^ \d) via F^; \d) \b) via T^i; \d) — |c) via F^q- 
Accordingly, using Eqs. (0) we get 



O'aa 


= -^LC^aa 


+ r_R0<^fe6 + ^RlO'cc 5 


(8a) 




= ^LCTaa " 


(XbO + ^in)crbb + ^RlO'dd , 


(8b) 


Occ 


= ^inO'bb - 


(Xl + ^Rl)(^cc + ^RO<^dd , 


(8c) 




= ^L(7cc - 


(XrO + ^Rl)o'dd ■ 


(8d) 



The total current / = Irq + Iri is given by Eq. (^, where Irq = /qTuq is the current fiowing 
to the right reservoir from the upper level, and Iri = fiTjn is the current fiowing from the 
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lower level (the "inelastic current"). Here /o = (Xfef, + add and /i = o"cc + o'dd are the electron 
densities on the upper and the lower levels, correspondingly. After solution of Eqs. (H) we 
obtain 

™ " Ft + aiVL + Tno) ' " + a(r^ + T^o) ' ^ ^ 

where Tr = + + ^in and a = riri„/[rRi(rL + Trq + Tri)]. 

This result is quite different from the stationary resonant current having been 

found within the one-electron framework 



^(one-el) ^ RO ^(one-el) ^L^in /'in"\ 

^RO — ' ^1 ~ p ■ ^ 

Notice that F^i does not enter Eq. ([l^), and therefore the stationary current in the one- 
electron approach does not depend on subsequent processes following the electron relaxation 
to the lower level. This follows from the irreversibility of the inelastic scattering, which 
disconnects the inelastic flux from the initial channel Hence, it is very remarkable 

that the F/ji-dependence of the resonant current appears in Eq. (^. It means that, when 
the exclusion principle is taken into account, the irreversibility does not totally disconnect 
the processes taking place before and after the relaxation. 

One can show |^ that the inelastic current Iri^~°^^ in the one-electron approach, Eq. (|lO|), 
can be represented as Iri'^~°^^ = fo^m, where /o is the electron density on the upper level. 
Therefore one may tempt to account for the exclusion principle by writing the inelastic 
current as Iri = /o(l — /i)rj„, where fi is the electron density on the lower level, Eq. @. 
However, one can easily check, by using Eq. (||), that this assumption does not hold. This 
example shows explicitly that ad hoc introduction of the additional factor (1 — /) cannot be 
correct in general for taking account of the Pauli principle in quantum transport PJ^,p!5|]. 

It is important to mention that the total current I = Irq + Eq. (^, always increases 
with Tin- It is similar to the result of the one-electron treatment, Eq. (p!0|). However, the 
account for the exclusion principle leads to restriction of the total resonant current / when 
Tin oo (as can be easily seen from Eq. (^). 



The most spectacular manifestation of the Pauh exclusion principle would take place in 
the resonant tunneling through double-well structures, in presence of inelastic transitions. 
As an example we consider the system shown in Fig. 2b, where a resonant current flows due 
to inelastic transitions from the upper to the lower levels in the left well. We investigate the 
case where the lower level in the left well, Ei, and the level E2 in the right well are aligned, i.e. 
El = E2, so the electrons may oscillate between the wells. As a result, the exclusion principle 
generates non-trivial quantum interference effects in the course of inelastic transitions from 
the upper level Eq to the lower one Ei 



Enumerating all possible electron states one can straightforwardly write down the system 
of linear coupled equations, Eqs. (H), for all transitions in the device. We do not present 
these equations here, since the procedure is exactly the same as in all the previous examples. 
Here we give only the final expression for the resonant current, Eq. (|^), when the inelastic 
rate Tin > Tl,r, f^: 

,^_Mr^r^ (11) 

Here the coefficients A4, Af and V are functions of the tunneling widths, Tl^h, and the 
coupling Q between the aligned levels, namely 

;w = ^]2(^2 +4^2) + 2^i^l^^ (i2a) 

Af = ^,(2^]2 + T,mT' - T,Tn) + ^ (r^ - ri + , (12b) 

V = TiTlTRin^ + r^) , (12c) 

where T = {Tl + Tr)/2 and Fi = {2Vl + Tr)/2. It follows from Eq. ( pT]) that, contrary to 
expectations, the resonant current decreases to zero when the inelastic rate Fj„ increases. 
Such a strange behavior of the tunneling (inelastic) current is a result of quantum interference 
effects due to the Pauli principle. Indeed, since the tunneling time from the left well to the 
right one is much larger than the time of inelastic transition from the upper level, an electron 
from the upper level should have arrived the second well "simultaneously" with the previous 
electron. It leads to an effective localization of electrons at the upper and the lower levels 
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in the left well, and therefore to the decrease of the inelastic resonant current. 

One can also consider a configuration where the level Ei in Fig. 2b is above the bottom 
of the Fermi sea in the left reservoir. Then the total resonant current has two components: 
the first one from the inelastic transitions from the upper level Eq, and the second one 
from the direct tunneling to the level Ei. In this case, too, the first component interferes 
destructively with the second one. As a result, the total resonant current increases when 
the first component is "switched off" (i.e. when Ep < Eq). 

All the discussed deviations from the single particle picture can be verified in experiments 
with semiconductor quantum wells. Yet, one should take into account that inelastic transi- 
tions with subsequent rescattering on impurities, with the momentum transfer to directions 
parallel to the layers, may decrease the expected deviations. 
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FIGURES 

FIG. 1. Resonant tunneling through (a) double-barrier and (b) triple-barrier heterostructures 
in absence of inelastic scattering. 

FIG. 2. Resonant tunneling through (a) double-barrier and (b) triple-barrier heterostructures 
in presence of inelastic scattering. 



12 



